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We show how in principle to construct analogs of quantum Hall edge states in "photonic crystals" 
made with non-reciprocal (Faraday-effect) media. These form "one-way waveguides" that allow 
electromagnetic energy to flow in one direction only. 
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In this letter, we describe a novel effect involving an 
interface between two magneto-optic photonic crystals 
(periodic "metamaterials" that transmit electromagnetic 
waves) which can theoretically act as a "one-way waveg- 
uide", i.e., a channel along which electromagnetic en- 
ergy can propagate in only a single direction, with no 
possibility of being back-scattered at bends or imper- 
fections. The unidirectional photonic modes confined 
to such interfaces are the direct analogs of the "chi- 
ral edge-states" of electrons in the quantum Hall effect 
(QHE) The key enabling ingredient is the presence 

of "non-reciprocal" (Faraday-effect) media that breaks 
time-reversal symmetry in the metamaterial. 

Just as in the electronic case, every two dimensional 
photonic band is characterized by a topological invariant 
known as the Chern number^, an integer that vanishes 
identically unless time-reversal symmetry is broken. If 
the material contains a photonic band gap (PBG), the 
Chern number, summed over all bands below the gap, 
plays a role similar to that of the same quantity summed 
over all occupied bands in the electronic case. In particu- 
lar, if the total Chern number changes across an interface 
separating two PBG media, there necessarily will occur 
states localized to the interface having a non-zero net 
current along the interfacep], Q. In the photonic case, 
such states would comprise our "one-way waveguide" . 

Such an interface between two PBG media can be re- 
alized as a domain wall in a 2D periodic photonic meta- 
material, across which the direction of the Faraday axis 
reverses. Unidirectional edge states are guaranteed in 
this system provided that the Faraday effect generates 
photonic bands with non-zero Chern numbers. Here, we 
construct photonic bands with non-zero Chern invariants 
in a hexagonal array of dielctric rods with a Faraday ef- 
fect. We then show that as a consequence of topology 
of the single-particle photon bands in the Brillouin zone, 
the edge states of light occur along domain walls (where 
the Faraday effect vanishes). 

It may seem surprising that the physics of the QHE 
can have analogs in photonic systems. The QHE is ex- 
hibited by incompressible quantum fluid states of elec- 
trons - conserved strongly-interacting charged fermions - 
in high magnetic fields, while photons are non-conserved 



neutral bosons which do not interact in linear media; fur- 
thermore, photonic bands can be described classically, in 
terms of Maxwell's equations. However, the integer QHE 
can in principle occur without any uniform magnetic flux 
density (just with broken time-reversal symmetry) as has 
explicitly shown by one of us in agraphene-like model 
of non-interacting Bloch electrons [6|]; thus Landau-level 
quantization is not an essential requirement for the quan- 
tum Hall effect. 

We have transcribed the key features of the electronic 
model of Ref. to the photonic context. The edge-states 
are a property of a one-particle eigenstate problem sim- 
ilar to the Maxwell normal-mode problem, so are repli- 
cated in the photonics problem. (The QHE itself has no 
photonic analog, as it follows from the Pauli principle of 
filling all one-particle states below the Fermi level.) 

The Maxwell normal-mode problem in loss-free linear 
media with spatially-periodic local frequency-dependent 
constitutive relations is a generalized self-consistent Hcr- 
mitian eigenproblem, somewhat different from the stan- 
dard Hermitian eigenproblem. The non-reciprocal parts 
of the local Hermitian permittivity and permeability ten- 
sors e{r,ui) and (i(r,u>) are odd imaginary functions of 
frequency, so frequency-dependence is unavoidable. The 
generalized eigenproblem has the structure 

C/t(fc)At/(fc)K(fc)> =«„(fc)B(w„(fc))|« n (fc)), (1) 

where U(k) is a unitary operator that defines the Bloch 
vector k; A and B(u>) are Hermitian operators, with 
the real-eigenvalue stability condition that the Hermi- 
tian operator Bq(lu) = (d/du))(ujB(v)) is positive defi- 
nite (this assumes that the periodic medium coupled to 
the electromagnetic fields has a linear response described 
by harmonic oscillator modes, none of which have natu- 
ral frequency u n (k) - a detailed derivation has been pre- 
sented in Ref. Q). The eigenfunctions (r\u n (k)} are the 
spatially-periodic factors of the Bloch states. The elec- 
tronic band-structure problem is a simplification of (fTj). 
with A replaced by the one-electron Hamiltonian, B by 
the identity operator 1, and w„ by the energy eigenvalue. 

In this formulation of Maxwell's equations, the eigen- 
function u n (k,r) = (r\u n (k)) is the 6-component vector 
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of complex electromagnetic fields 



u n (k,r) 



E n (k,r) 
H n (fc,r) 



(2) 



In this basis, U(k, r) — expik-r and A = —iJ a V a (with 
V a = d/dr a , and repeated indices summed), where 
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iL a 
-iL a 
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e(r, uj) 







M(r,w) 



(3) 



here L a are the 3x3 spin-1 matrices, (L°) ac = ie . If the 
physical electromagnetic fields are given by the real parts 
of u n (k, r) expi(k ■ r — u> n (k)t), the spatially-periodic 
time-averaged energy density and energy current are the 
quadratic forms (u* , Bo(u> n )u n ), (u*, J a u n ). For our 
purposes, the key photonic band-structure quantity is the 
Berry connection A„(k), a real function of k given by 



AO, 



t |B K)|V> K ) - (Vgu n |Bo(^n)|Mn) 

2i(u n \B a (uj n )\u n ) 



(4) 



where = d/dk a is the fc-space derivative. We ob- 
tained ([?]) as a generalization of the B = 1 expression[5( 
by deriving |1| from a standard Hermitian eigenprob- 
lem where electromagnetic fields are explicitly coupled 
to harmonic-oscillator degrees of freedom of the medium 

The solution of the normal-mode eigenproblem 
only determines u n (k,r) up to an arbitrary k- 
dependent phase factor; if the replacement u n (k, r) — > 
u n (k,r)expixn(k) is made, A„{k) -> _4£(fc) + V£xn( fc )- 
The Berry connection is a "gauge-dependent" analog 
of the electromagnetic vector potential; the associated 
gauge-invariant function (analogous to the magnetic flux 
density) is the fc-space Berry curvature J-„ b (k) — V^,A b n — 
\7 b k An- The Berry phase expi</>„(r) = expi^A^dka 
associated with adiabatic evolution around a closed 
path r (here in fc-space) is the gauge-invariant analog of 
the Bohm-Aharonov phase factor, and can be expressed 
in terms of the integral of over a surface bounded by 

ri- 

The Berry curvature satisfies a fc-space analog of the 
Gauss law, except that "monopole" singularities emitting 
total "Berry flux" ±2tt can occur at fc-space points where 
there are "accidental degeneracies" between bands (this 
quantization of the monopole charge ensures that the ex- 
pression for the Berry phase in terms of Berry curvature 
on a surface bounded by T is independent of how that 
surface is chosen^). The integer Chern invariant associ- 
ated with any compact surface (2-manifold) £ in fc-space 
is 
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dk a A dkb T', 



ah 



(5) 



In the case of a 2D band-structure, S may be taken to be 
the 2D Brillouin zone (BZ) itself, and eft is a property 



of the 2D bandQ. If time-reversal symmetry is unbroken, 
F% (— k) = —J^ (k), and Chern numbers vanish. 

We now wish to construct a 2D photonic bandstruc- 
ture where some bands have a non-zero Chern number. 
The key idea is to start with a bandstructure that has 
both time-reversal symmetry and inversion symmetry, 
which allows the existence of pairs of "Dirac points" in 
the 2D BZ. These are isolated points where two bands 
become degenerate, but split apart with a linear disper- 
sion (resembling that of the massless Dirac equation) for 
nearby Bloch vectors. Dirac points are generically al- 
lowed because if inversion symmetry is present, !F^ b (~k) 
= J r ^ b (k); in combination with time-reversal symmetry, 
this means that JF° b (fc) = 0. It is then possible to chose 
a phase convention such that the eigenfunctions u n (k, r) 
of ((T|) are real for all k. The eigenproblem is then of the 
real-symmetric type, where it is possible to find an "ac- 
cidental" degeneracy between two bands by varying just 
two parameters (in this case, the 2D Bloch vector k); 
in contrast, in the general complex-Hermitian case, three 
parameters must be varied to find a degeneracy, which 
cannot be done by merely "fine-tuning" a 2D k. 

Dirac points can exist in a 2D bandstructure with both 
spatial-inversion and time-reversal symmetry, but a gap 
opens if either symmetry is broken. While breaking of 
inversion symmetry leads to non-zero Berry curvature 
(and hence corrections to the "semiclassical" equations 
for the trajectories of light rays in adiabatically-varying 
media|8|), it does not lead to non-trivial topology of the 
bands. In contrast, when a gap opens at Dirac points due 
to time-reversal breaking, the two bands that split apart 
inevitably acquire non-zero Chern numbers. 

We can now give an in-principle demonstration that 
"one-way waveguides" can be constructed using non- 
reciprocal photonic crystals. Consider a system with 
a uniform isotropic permeability tensor /io<5 ab , and 
an isotropic but spatially-varying permittivity tensor 
e(r)e a 6 ab , with 



e(r) = e(l + AVfc(r)), V G {r) 



2^cos(G„-r), (6) 



where G n , n = 1,2,3, are three equal-length reciprocal 
vectors in the xy plane, rotated 120° relative to each 
other. For small A, this problem can be solved ana- 
lytically in a "nearly-free-photon" approximation. This 
system has continuous translational invariance in the z- 
direction, and we will restrict attention to wavenumbers 
k z = 0, with Bloch vector k in the xy plane. The electro- 
magnetic fields then separate into decoupled "TE" and 
"TM" sets, {E x ,E y ,H z } and {H x ,H y ,E z }; we special- 
ize to the TE set. The six corners of the (first) BZ arc 
at ±K n , where K\ — (Gi — Gs)/3, etc., and \K\ = 
\G\/y/3; since K2 — K\ = G3, the three wavevectors 
K.i are equivalent as Bloch vectors. To leading order 
in A and the deviation 5k = k — Ki of the 2D Bloch 
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vector from the the BZ corner, the three "free photon" 
TE plane waves with speed cq split into a "Dirac-point" 
doublet with lu = lu d ± vo\Sk\ + 0(\6k\ 2 ), where lu d = 
coli^Kl - A/4 + (9(A) 2 ), v D = c o /2 + 0(A), and a singlet 
lu = lu q + 0(\Sk\ 2 ), lu = c \K\(l + A/2 + <3(A 2 )). 

We now perturb the Dirac points by a Faraday term 
(which explicitly breaks time-reversal symmetry), with 
an axis normal to the xy plane, added to the permittivity 
tensor: e xy = —e yx = ieoerj(r,ui), where 

V(r,uj) = 770 M +r}x(u))V G {r); (7) 

770 (w), Tj\{u>) are real odd functions of lu. We assume that, 
for lo sa lud, \r]o(uj)\,\r]i((j)\ <C |A| <C 1, with negligible 
frequency-dependence. The Dirac points now split, with 
dispersion lu — lud ±V]j(\5k\ 2 + n 2 ) 1 / 2 , where, to leading 
order in -q, k = \K\(^T)i(oj d ) - 3A?7o(a; D )). 

For small k, the Berry curvatures of the upper and 
lower k z = bands near the split Dirac points are 

FZ v {5k) = ±\n(\8k\ 2 + k 2 Y V2 . (8) 

There is a total integrated Berry curvature of ±7r near 
each Dirac point, giving total Chern numbers ±1 for the 
split bands. By inversion symmetry, the Berry curvatures 
at the two Dirac points have the same sign; if the gap 
was opened by broken inversion symmetry, with unbro- 
ken time-reversal invariance, they would have opposite 
sign, and the Chern number would vanish. 

We now consider an adiabatically spatially-varying 
Faraday term parameterized by a n(r) that is positive 
in some regions and negative in other regions. The split- 
ting of the Dirac points vanishes locally on the line where 
n(r) = 0. It is necessary that, in the perfectly periodic 
structure with k = 0, there arc no photonic modes at 
other Bloch vectors that are degenerate with the modes 
at the Dirac points. 

Such frequency- isolation of the Dirac points cannot oc- 
cur in the weak-coupling "nearly- free photon" limit, but 
can be achieved, at least for k z = modes, in hexagonal 
arrays of infinitely long dielectric rods parallel to the z 
axis. An example can be seen in Fig. (la) of Ref. [sj. That 
figure was exhibited to demonstrate a frequency gap be- 
tween the first and second TE bands, but incidentally 
also shows that the second and third TE bands are sep- 
arated by a substantial gap except in the vicinity of the 
BZ corners, where they touch at Dirac points. The cor- 
responding TM bands were not given in Ref. [9], but we 
found that the Dirac-point frequency u>d is also inside 
a large gap of the TM spectrum (see Fig.JT])). When a 
Faraday term is added, the bands forming the Dirac point 
in Fig. ([I]) split apart, and each now non-degenerate band 
will have associated with it a non-zero Chern number (see 
Ref. (3). 

The Faraday effect incorporated to the hexagonal ar- 
ray of rods explicitly breaks time-reversal symmetry on 
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FIG. 1: Photon bands for k z — electromagnetic waves prop- 
agating normal to the axis of a hexagonal 2D array of cylin- 
drical dielectric rods; a is the lattice constant. As in Fig. (la) 
of Ref. 0|, the rods fill a fraction / = 0.431 of the volume, 
with dielectric constant e = 14, and are embedded in an e = 
1 background. The lowest five 2D bands are well-separated 
from higher bands, except near a pair of "Dirac points" at 
the two distinct Brillouin zone corners (J). 



the scale of the unit cell of the metamaterial: the permit- 
tivity tensor acquires an imaginary, off-diagonal compo- 
nent having the periodicity of the unit cell, as described 
above. A hexagonal array consisting of a material hav- 
ing a large Verdet coefficient, such as a rare-earth garnet 
with ferromagnetically-ordered domains would give rise 
to such an effect. 

While these k z = Dirac-point modes are not degen- 
erate with any other k z = modes, they are degenerate 
with k z modes. To fully achieve a "one-way" edge- 
mode structure, the light must also be confined in the 
z-direction, with Dirac points at a frequency that is non- 
degenerate with any other modes. To design such struc- 
tures, it will be necessary to vary the filling factor of the 
rods along the z-direction, so that light remains confined 
to regions of relatively larger filling factors. The technical 
challenge would be to vary the filling factors without in- 
troducing any modes into the bulk TE gaps surrounding 
the Dirac points. 

Let \u a (±kr>)), & — ±, be the degenerate solutions 
of ([TJ at a pair of isolated Dirac points, normalized 
so (u cr (±k D )\B (uj D )\u rT/ (±k D )) = B Q S aa >. Now add 
a Faraday perturbation SB(t,lu): in degenerate pertur- 
bation theory, normal modes with small 6lu — lu — lud 
have the form ^ a ± ipf(r)U(±kri, r)u a (±kn 7 r). For 
slow spatial variation, there is negligible mixing between 
modes at different Dirac points, and ip^(r) is the solution 
of 

Y, (-iJlVa - u D SB(r))± a , ^±(r) - SuB Q ^(r), 

(9) 
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where JJ and SB(r) are 2x2 matrices given by 

(Jl)± a , = (u a (±k D )\J a \u a ,(±k D )), a = x,y, 
( 5B ( r ))**' = (ua{±k D )\5B{r,oj D )\u a ,{±k D )).{lQ) 

For a straight-line interface, this equation has the form 
VjjK\ip) = 8u)\ip), with Vjj > 0, and 

K = -icr x V x + Sk ll( T y + k(x)(t z , (11) 

where er° are Pauli matrices. Here ko y + 8k\\ is the con- 
served Bloch vector parallel to the interface; we take k{x) 
to be monotonic, with k(x) — > ±k°° as x — > ±oo. 

It is instructive to first consider the exactly-solvable 
case k(x) — k°° t&nh(x / , £ > 0, where K 2 is essen- 
tially the integrable Poschl- Teller Hamiltonian[10]. The 
spectrum of modes bound to the interface is 

u)o(6k\\) = lod + s K v D 8k\\, s K ee sgn(ft°°), (12a) 

wn±(<%) = ijj D ±v D \8k\ + K i n j , n>0, (12b) 

with | K n \ < | k°° | ; for the integrable model, k 2 is given 
by 2n|K°°|/^, n < |«°°|£/2. There is always a unidirec- 
tional n — mode with speed vd and a direction deter- 
mined by the sign of k°°; in the small-£ (or sharp- wall) 
limit < 2, this is the only interface mode. 

Let 4>(k 2 ) be the dimensionless cirCci in the x-k% 
phase-plane enclosed by a closed constant-frequency orbit 
(k x ) 2 + ( K {x)) 2 — k 2 < |k°°| 2 , corresponding to a bound 
state. For the integrable model, this has the simple form 
</>(k 2 ) = ttk 2 ^/\k°°\; the n > bidirectional modes thus 
satisfy a constructive-interference condition 

4>(kI) = 27m. (13) 

This contrasts with the usual "semiclassical" condition <p 
= 2it{n + |); the change is needed for the n = "zero 
mode" (|12al) to exist, and can be interpreted as deriving 
from an extra Berry phase factor of —1 because the orbit 
encloses a Dirac degeneracy point at (x, k x ) = (0, 0). For 
general k(x), the n = eigenfunction is 

i/j°(r) oc v3 CT (s K )exp (i6k\\y - s K J K(x')dx^j . (14) 

er v ip(s) = s<p(s). For slowly- varying k(x), the condition 
(|13| will determine k 2 for any n > interface modes. 

Since there are two Dirac points, there are two such 
unidirectional edge modes at a boundary across which 
the Faraday axis reverses. The crucial feature is that 
both modes propagate in the same direction, and cannot 
disappear, even if the interface becomes sharp, bent, or 
disordered. As in the QHE, the difference between the 
number of modes moving in the two directions along the 
interface is topologically determined by the difference of 
the total Chern number of bands at frequencies below 



the bulk photonic band gap in the regions on either side 
of the interface; in this case lAC^'l = 2. 

For \8ui\ < a Faraday interface has no counter- 

propagating modes into which elastic backscattering can 
take place, so the "one-way waveguide" that it forms is 
immune to localization effects, just like electronic trans- 
port in the QHE. In the QHE, the number of electrons is 
strictly conserved; in photonics, the photons only propa- 
gate ballistically if absorption and non-linear effects are 
absent. These effects do allow degradation of the elec- 
tromagnetic energy current flowing along the interface, 
so the analogy with the QHE is not perfect. 

Even if a 2D metamaterial with isolated Dirac points 
can be designed, the problem of finding a suitable 
magneto-optic material to provide the Faraday effect 
must be addressed. The effect must be large enough to 
induce a gap that overcomes the effect of any inversion- 
symmetry breaking. The parameter is the inverse 
length that controls the width of the unidirectionally- 
propagating channel (and the unidirectional frequency 
range); in order to keep the wave confined to the inter- 
face, and prevent leakage, the Faraday coupling must be 
strong enough so that this width is significantly smaller 
than the physical dimensions of the sample of metama- 
terial. 

In summary, we have shown that analogs of quantum 
Hall effect edge modes can in principle occur in two di- 
mensional photonic crystals with broken time-reversal 
symmetry. The electromagnetic energy in these modes 
travel in a single direction. Explicit theoretical examples 
of such modes have been constructed in Ref. O]). Such 
quasi-lossless unidirectional channels are a novel possibil- 
ity that might one day be physically realized in "photonic 
metamaterials" with non- reciprocal constituents. 
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